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We describe a measurement of the direct CP asymmetry between inclusive 6 sj and b ^ 
decays. This asymmetry is expected to be less than 0.01 in the Standard Model, but could be 
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enhanced up to about 0.10 by new physics contributions. We use a sample of 89 million BB 
pairs recorded with the BABAR detector at PEP-II, from which we reconstruct a set of 12 exclusive 
b ^ sj final states containing one charged or neutral kaon and one to three pious. We measure an 
asymmetry of Acp{b s'y) = 0.025 ± 0.050{stat.) ± 0.015{syst.), corresponding to an allowed range 
of -0.06 < Acp{b S'y) < +0.11 at 90% confidence level. 

PACS numbers: 13.35.Dx, 14.60.Fg, 11.30.Hv 



The inclusive decay 6 — > S7 is a flavor-changing neu- 
tral current process described by a radiative penguin 
loop diagram. The world average branching fraction is 
(3.5 ± 0.5) X 10"'' 1] in good agreement with recent the- 
oretical predictions (2]. Earlier experimental values of 
the branching fraction have been used to constrain new 
physics beyond the Standard Model 3]. A measurement 
of the direct CP asymmetry between 6 — > 57 and b —^sj 
decays provides an independent and significant test of 
these predictions. In the Standard Model the dominant 
loop contribution contains a top quark, with other contri- 
butions being suppressed by CKM factors and the GIM 
mechanism. The lack of interference between comparable 
amplitude contributions leads to a rather small predicted 
asymmetry 4]: 



^CP 



T{b — > 57) — r(6 — ^ 57) 
T{b 57) + T(b ^ 57) 



OmAAt^oT, (1) 



which has little sensitivity to the photon energy cut-off or 
to the distribution of hadronic final states. The dominant 
errors are due to the uncertainty of the charm quark mass 
and the choice of the perturbative scale. The inclusion of 
contributions to the loop beyond the Standard Model can 
increase the predicted asymmetry up to about 0.10 0. 

There is a previous measurement of direct CP asym- 
metry in a sum of 6 ^ 57 and b ^ d'y decays. In 
the Standard Model, the total of the b ^ sj and b ^ d'y 
asymmetries is exactly zero in the U-spin symmetry limit, 
rud = ms, as a consequence of CKM unitarity The 
measurement in Ref. p gives —0.27 < 0.965 x Acp{b 
s-f) + 0.02 X Acrib^ d-y) <0 AO. _ 

We use a sample of (88.9±1.0) x 10^ BB pairs collected 
at the T(4S) resonance with the BABAR detector at the 
PEP-II asymmetric e^e~ collider. A detailed description 
of the detector can be found elsewhere . For this anal- 
ysis the most important detector elements are the forty- 
layer drift chamber, situated in a 1.5 T solenoidal mag- 
netic field, which measures charged particle momenta, 
the CsI(Tl) electromagnetic calorimeter, which measures 
the energies of the photons, and the detector of inter- 
nally refiected Cherenkov light (DIRC), which is used to 
identify charged kaons. 

We reconstruct b ^ sj decays as the sum of twelve 
exclusive final states: 

B^ —f K^TT^J, K^TT^TT^J, K^Tt'^Tt'^J, K^TT^TT^Tt'^J 

_R'"7r+7, _ftr"7r+7r'^7, iir~7r~'"7r'^7r*'7, _Rr~7r+7r"7r+7 
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and measure the yield asymmetry with respect to their 
charge conjugate decays b 57. The identification of 
charged kaons removes b ^ d'y decays. We do not use 
decays to final states with iiTg to determine the direct 
CP asymmetry, since these are not fiavor-specific, but 
we study them to understand systematic effects. 

The high energy photon is detected from an isolated 
energy cluster in the calorimeter, with shape consistent 
with a single photon, and energy E* > 1.8 GeV in the 
e+e~ center-of-mass frame. A veto is applied to the high 
energy photons that combined with another photon form 
either a 7r° within the mass range 117-150 MeV/c^ or an 
77 within the mass range 524-566 MeV/c^. 

Neutral kaons are reconstructed as K*^ tt+tt" can- 
didates with an invariant mass within 9 MeV/c^ of the 
nominal mass , and a transverse flight distance > 2 mm 
from the primary event vertex. Charged kaons are tracks 
identified as kaons from information in the DIRC. The re- 
maining tracks are considered to be charged pions. Both 
charged and neutral kaons are required to have a labora- 
tory momentum > 0.7 GeV/c. Above this threshold the 
rate for charged pions to be mis-identified as kaons is < 
2.0%. 

Neutral pions are reconstructed from pairs of photons 
with energies > 30 MeV. A tt'^ mass cut is applied be- 
tween 117 and 150 MeV/c^. Charged and neutral pions 
are required to have laboratory momenta > 0.5, 0.3 or 
0.2 GeV/c for states with 1, 2 or 3 pions, respectively, to 
reject combinatoric background. 

The mass of the hadronic system, X^, formed from the 
kaon and pions is required to be between 0.6 GeV/c^ and 
2.3 GeV/c^, corresponding to a photon energy threshold 
> 2.14 GeV in the B rest frame. 

The signal Monte Carlo sample is generated according 
to Ref. y, which predicts that (83 ± 5)% of the b 
S7 spectrum is above our photon energy threshold. We 
use JETSET to hadronize the system of the strange 
and spectator quarks. Within the selected hadronic mass 
range, the twelve final states constitute 48% of the total 
rate. If we also include the 5° decays to and equate 
the decays to with those to this increases to 73% 
of the total rate. As a part of our analysis, we check the 
dependence of the asymmetry on the hadronic mass and 
final state. 

Most of the background in this analysis arises from 
continuum production of a high energy photon, either by 
initial state radiation, or from the decays of tt*' and r] 
mesons. We remove 86% of these backgrounds by selec- 
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tions on the angle between the thrust axis of the B meson 
candidate and the thrust axis of all the other particles of 
the event, |cos0^| < 0.80, and the angle between the 
B candidate and the beam axis, |cos0^| < 0.80, both 
defined in the e+e" center-of-mass system. We then 
use a neural network to combine information from a set 
of event shape variables, including a set of energy flow 
cones. This halves the continuum background compared 
to our initial selection. 

In 12% of the signal events, we can identify an electron 
or muon from the decay of the other B [lo| . This is a very 
effective signature for removing continuum background, 
so the remaining background in this sample comes mostly 
from other B decays. We present separately our results 
for the sample of events which are lepton-tagged. 

Exclusive & — > S7 decays are characterized by two 
kinematic variables: the beam-energy substituted mass, 
niES ~ vTnT^T^P'^-P^, and the energy difference be- 
tween the B candidate and the beam energy, Ai? — 
Eg — where E*^ and p*g are the energy and mo- 

mentum of the B candidate in the e+e~ center-of-mass 
frame, and ^/s is the total center-of-mass energy. We 
require candidates to have |Ai5| < 0.10 GeV, and re- 
move multiple candidates in each event by selecting the 
one with the smallest value of |Ai?|. This technique is 
> 90% efficient when the true 6 — > 57 decay is among 
the reconstructed candidates. We then fit the niES dis- 
tribution between 5.22 and 5.29 GeV/c^ to extract the 
signal yield. When calculating tobs, the value oi p*^ is 
corrected for the tail of the high energy photon response 
function by scaling the measured E* to the value that 
would give Ai? — 0, the value expected for true signal. 

In order to fit the toes distribution in data, we need 
to understand the different components of the signal and 
background events. We have identified the following four 
contributions as shown in Figure^ The signal events are 
described by a Crystal Ball function with a resolu- 
tion a{mEs) = 2.2 MeV/c^. The continuum background 
is described by an ARGUS shape 0|, which is cross- 
checked by a fit to a sample of 9.6/6~^ of data taken 
40 MeV/c^ below the r(4S) resonance. We use a B'B 
Monte Carlo sample to model the background from B 
decays other than b — + S7, which is significant for Xs 
masses above 1.9 GeV/c^. This background is described 
by the sum of an ARGUS shape and a peaking compo- 
nent which is modelled by the signal shape. 

The last background component is cross-feed from in- 
correctly reconstructed 6 — > 57 events. This is modelled 
by the signal Monte Carlo sample, where we identify 
events reconstructed in the wrong final state. Cross-feed 
occurs when the true 6 — > 57 decay is not among the re- 
constructed candidates, or in a multiple candidate event 
when the wrong candidate is chosen. The shape of the 
cross-feed is described by the sum of an ARGUS shape 
and a peaking signal shape. We regard cross-feed as a 
background to be subtracted. 
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FIG. 1: Monte Carlo simulations of the four contributions 
to the beam-energy substituted mass distribution of events 
selected as & ^ 57, with the corresponding fits: (a) signal, 
(b) continuum, (c) BB decays and (d) cross-feed. The plots 
are normalized to the luminosity of our data sample. 



We fit the data toes distributions separately for 
each fiavor. For the total sample, the fit function is 
parametrized by two ARGUS shapes and a Crystal Ball 
function. One ARGUS shape is fixed to be as the con- 
tinuum ARGUS shape, while the other one is free to rep- 
resent the sum of the non-peaking BB and cross-feed 
backgrounds. The Crystal Ball function fits the combi- 
nation of the peaking components. For the lepton-tagged 
sample, we use only one free ARGUS shape and a Crystal 
Ball function. In all cases we use an unbinned maximum 
likelihood fit. The fitting technique has been validated 
with a large sample of Monte Carlo simulated events. 
In Figure El we present the final fits to the to ^5 distri- 
butions for 5 — + S7 and 6 — > 57 events. The lower plots 
are for the lepton-tagged sample. All the fits have per 
degree-of-freedom close to 1 , if we make a fit to a binned 
distribution as shown in Figure [21 The sum of events in 
the h and b peaks is 1644 ± 72, of which 201 ± 18 are 
lepton-tagged. To get the true signal yields these have 
to be corrected for the predicted yield of peaking BB 
and cross-feed backgrounds from Monte Carlo samples 
(see Figure^, which is 88± 27, where 10 ± 8 are lepton- 
tagged. 

The direct CP asymmetry is calculated from: 



A 



CP 



1 

W) 



{n — n) 
[n + n) 



aDET 
^CP 



(2) 



where n and n are the numbers of observed b ^ sj 
and 6 ^ 57 events after the peaking background is 
subtracted, AD = 2{w — w) is the difference in the 
wrong flavor-fraction between b and b decays, and (D) = 
1 — {w + ad) is the dilution factor from the average 
wrong flavor-fraction. A^p"^ is the flavor-asymmetry 
of the detector. We find AD = 0.001 ± 0.002 and 
(D) = 0.989 ± 0.001 from Monte Carlo samples. The 
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FIG. 2: Fits to the beam-energy substituted mass distribu- 
tions in data events for: (a) all b —> S7, (b) all b — > 57, (c) 
lepton-tagged b —> and (d) lepton-tagged 6 ^ 37 decays. 
Contributions are shown from peaking Crystal Ball (dotted- 
dashed), fixed continuum ARGUS shape (dotted) and free 
BB and cross-feed ARGUS shape (dashed). 

small wrong-flavor fraction is due to charged pions mis- 
identified as charged kaons. 

We need to correct the measured value of Acp for the 
flavor-asymmetry of the detector A^p^ . While it is 
known that the kaon-nucleon cross-sections are asym- 
metric at low momenta, there are few accurate mea- 
surements Q. This means that our Monte Carlo sam- 
ple is not expected to model correctly the asymmetries 
due to the interactions of kaons with the inner part of 
the detector. The kaon identification efficiency of the 
DIRC for reconstructed tracks is measured with a con- 
trol sample of kaons from D* decays. Averaging over 
the kaon spectrum in 6 — *■ 57 events we obtain a small 
asymmetry of —0.002 ± 0.001 from particle identifica- 
tion. We measure the overall detector asymmetry of the 
data events in our wies and H^E sidebands, increasing 
the statistics by removing the neural network cut. Most 
of these events are from the continuum, where we do 
not expect any physics mechanism to generate a flavor- 
asymmetry. We observe a significant asymmetry for kaon 
momenta below 1 GeV/ c. The asymmetry as a function 
of the kaon momentum is applied to the signal Monte 
Carlo to determine what shift should be applied to the 
data. This gives an overall flavor-asymmetry correction 
A^^'^ = -0.014 ±0.015. 

Table ^ presents the measured signal yields and cor- 
rected CP asymmetries. The lepton-tagged results are 
consistent with the results for the total sample. We di- 
vide the total sample into four bins in Xs mass, and 
observe no significant mass dependence of the asymme- 
try. The first bin corresponds to the if* (892) resonance, 
for which the world average asymmetry from studies 
of exclusive B K*^ decays is Acp{B — > K*^) — 
—0.005 ± 0.037 P|. Our result is consistent with this 



average. 

We divide our total sample into three types of decay 
mode: B°(B°) , B^ and B± ^ K^. We 

observe a discrepancy of 2.3cr between the two B^ cate- 
gories which we regard as a statistical fluctuation, since 
it is not correlated with a specific final state or hadronic 
mass bin. The combination of the B^ samples is consis- 
tent with a null asymmetry, as is the B'^ sample. 

TABLE 1: Signal yields and CP asymmetries for total and 
only lepton-tagged event samples. The total sample is also 
divided up into four bins in Xs mass in GeV/c^ , and into three 
types of decay modes. The errors on n and n are statistical 
only, while for Acp we quote the additional systematic error 
from the detector asymmetry. 



Sample 


n 


n 


Acp 


Total Sample 


787 ± 54 


769 ± 54 


0.025 ± 0.050 ± 0.015 


Lepton-tagged 


91 ± 14 


100 ± 13 


-0.04 ± 0.10 ± 0.02 


Afxs=0.6-l.l 
Mxs = l. 1-1.5 
Mxs=1.5-1.9 
Mxs=1.9-2.3 


378 ± 32 
162 ± 22 
139 ± 19 
101 ± 29 


396 ± 33 
136 ± 23 
124 ± 21 
67 ± 36 


0.003 ± 0.059 ± 0.015 
0.11 ± 0.11 ± 0.02 
0.07 ± 0.11 ± 0.03 
0.23 ± 0.30 ± 0.04 


B^^K^ 


455 ± 36 
229 ± 31 
100 ± 24 


447 ± 38 
148 ± 30 
166 ± 25 


0.015 ± 0.059 ± 0.014 
0.22 ± 0.12 ± 0.02 
-0.20 ± 0.14 ± 0.03 



The dominant systematic error in our measurement is 
the uncertainty of 0.015 in the fiavor-asymmetry of the 
detector. For the lepton-tagged sample we add an ad- 
ditional systematic uncertainty of 0.010 to account for 
a possible charge asymmetry in the lepton tagging ef- 
ficiericy. This is derived from studies of control sam- 
ples [13. 

We have tested the effect of possible flavor asymme- 
tries in the peaking cross-feed and BB backgrounds by 
varying them within the current experimental bounds 
(90% C.L.). We added a 0.10 asymmetry to the cross- 
feed events, and a 0.02 asymmetry to the peaking back- 
ground from BB decays, which comes primarily from 
B D^*^p decays. The change in our measured asym- 
metry due to these changes in the cross-feed and BB 
flavor-asymmetries is 0.004, which gives a negligible con- 
tribution to the error. 

We have checked that the parameters of the ARGUS 
shapes and Crystal Ball functions are the same for both 
flavors within Icr, so the detector asymmetry is simply an 
overall normalization difference between the two samples. 
We have also checked that the neural net distributions for 
signal and continuum background are flavor-symmetric. 

Our estimates of the cross-feed background and the de- 
tector asymmetry correction, Aqp'^ , depend on the mix 
of final states in our signal Monte Carlo sample. We 
check these, also using information from B'^ decays to 
final states with K^, by varying the ratios of final states 
with or if^, and tt" to tt"*" measured in our data 
by ±3(7. Note that the measured ratios are consistent 
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with our signal Monte Carlo. Changing the ratios has no 
significant effect on the cross-feed or the detector asym- 
metry correction. 

Our final result for the direct CP asymmetry in 6 — > 57 
is Acp = 0.025 ± 0.050 ± 0.015 for the total sample, 
and Acp = -0.04 ± 0.10 ± 0.02 for the lepton-tagged 
sample. The total sample provides the best constraint, 
-0.06 < Acp < +0.11 at 90% confidence level. 
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